The gy a7 - coupling constants in hadron 
electrodynamics 
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Abstract 



> 

Recent measurements of the branching ratios for the decays to — ► 7r°7r°7 and 
p° — > 7r°7r°7 lead to coupling constants of the Vo"7-interaction (1/ is a vector meson) 
one order of magnitude smaller than previously assumed to describe the threshold 
^ cross section of 7 + p — > p + p° and the HERMES effect. The new gy al couplings 

are in contradiction with the predictions of the QCD sum rules, but are in good 
agreement with VDM-estimation of the a — ► 77-width. 



1 Introduction 



The coupling constants gva-y, V = p or u are important ingredients for the 
theoretical analysis of many different hadronic electromagnetic processes. If 
m a < m v , the decays V — > a + 7 can occur directly- with radiation of electric 
dipole photons. In this respect, these decays are essentially different from the 
decays V — > P + 7, P = tt, r], rj' -with radiation of magnetic dipole photons. 
In the quark models, these last decays are induced by the quark magnetic 
moment, with transition S — 1 — > S — 0, where S is the total spin of the 
gg-system (in the corresponding meson), whereas the decays V — > a + 7 are 
induced by the internal motion of quarks, with transition S — 1 — ► S — 1 [1]. 

Let us mention the main applications of the ^y^-coupling constants: 
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• The a-exchange for the vector meson production [2-6] j + p — > p + p°, near 
threshold. 

• The electromagnetic transition 7* — > a + cj, in the space-like region of the 
virtual photon 7*four-momenta, enters in the calculation of meson exchange 
currents, in particular in the analysis of the deuteron electromagnetic form 
factors at large momentum transfer [7-9]. 

• The branching ratio of rare radiative decays of vector mesons as p° — > 7r°7r°7 
and uj — > 7r°7r°7 are controlled by the gvv 7 -coupling constants [10,11]. 

• The HERMES effect [12], concerning the inclusive electroproduction cross 
section on light nuclei can be explained under specific assumptions [13] 
about the absolute value of the g pal and g^ al coupling constants and the 
corresponding electromagnetic form factors. 

• The coupling constant of the a — > 77-vertex, which can be calculated on 
the basis of the gy ai couplings, is important for the analysis of real and 
virtual Compton scattering on nucleons [14-16]. 

• The exact value of the <?y CT7 coupling constant, in principle, allows one to 
constrain the g^AW-coupling, which describes the scalar exchange in the 
^-potential [17]. 

• Finally, the gy G1 coupling constants have an implicit theoretical interest, 
being fundamental coupling constants of hadron electrodynamics. 

Attempts to estimate these coupling constants from QC-D-sum rules, give rel- 
atively large absolute values [18]. Phenomenological methods based on the 
existing experimental data about different electromagnetic processes need se- 
rious additional assumptions. For example, in the framework of cr-exchange 
for 7 + p — > p + p°, it is possible to estimate the product g pa ^gaNN, assuming 
a definite form of the phenomenological form factors, which insure the correct 
behavior of the differential cross section da/dt, in the near-threshold region, 
and depend on ad-hoc cut-off parameters. Taking g% NN /4ir = 8, one obtains 
g pal = 2.71 [2]. Considering only cr-exchange in 7 + p — > p + p°, the sign of 
this constant can not be determined, but some polarization observables are 
sensitive to this sign [19]. Other possible contributions to the matrix element 
for the process 7 + p — > p + p°, such as iV*-excitation, iV-exchange in s- and 
li-channel are neglected in such oversimplified consideration [2]. 

Also the estimation of the discussed coupling constants done in framework of 
the HERMES effect [13] can not be considered as direct and model indepen- 
dent. 

Another source of information about the p<T7-vertex is the radiative decay 
p° — > 7r + 7r~7-with relatively large branching ratio [20]. Namely this decay has 
been considered to favor large values of g pay in comparison with 
of the corresponding widths T(u — > 7r + 7r~7) Y(p — > 7r + 7r~7) [2]. However it 
is necessary to point out that for the decay p° — > 7r + 7r _ 7 the bremsstrahlung 
mechanism dominates, therefore the a-contribution, being of the order of the 
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error bars, can not give the g pai with good accuracy. Note also that the precise 
determination of g pal depends essentially on the mass and the width of the 
cr-meson. In [11] different values of g pal could been found, with the general 
conclusion that the g paJ value from p° — > 7r + 7r _ 7 is not in contradiction with 
the estimation done in [2]. 

Finally, having a definite value of g pai , it is possible to estimate the branching 
ratio for p — > 7r°7r°7, in framework of the effective Lagrangian approach. In 
this case, the absence of the bremsstrahlung mechanism results in the fact 
that the a-exchange mechanism is important. But the predicted BR(p° — > 
7r°7r°7) ~ (45 — 200) • 10~ 5 , with g pay = 3 5, were too large in comparison 
with other theoretical expectations for this decay [21]. 

Moreover, and this will be the important background of this paper, large values 
of 

9paj contradict essentially the recent experimental results from Novosibirsk, 
concerning the direct measurement of the decay p° — > 7r°7r°7 [22]: 

Br(p -> 7r°7r°7) = (4.1±J;g ± 0.3) • lO" 5 , 

Br(p -> (77 -> 7T°A) = (1.9±g;l ± 0.4) • 10~ 5 . 

The main aim of this paper is to estimate the coupling constant g pai on the 
basis of these new important experimental data, which is, in our opinion, the 
most straightforward way. 



2 The (j-contribution to the p — > 7r° + 7r° + 7-decay 

A crude estimation of the g pal coupling constant can be easily obtained from 
the experimental data about p° — > (X7 — > 7r°7r°7 [22], under the assumption 
that the cr-mass is smaller than the p- meson mass, so that the decay p — > c + 7 
is allowed. 

Neglecting for a moment the a-width (which is an evident oversimplification 
of reality, as IV = (0.6 -7- 1) GeV [20]), we can find: 

r(p°-(77-7rV7) = ^r( A> -(77), (1) 

taking into account the identity of the neutral 7r°-mesons, produced in p° — > 
7r°7r°7, and the isotopic relation: g{a — > 7r°7r°) = g(a — > 7r + 7r _ ). 

The matrix element of the decay V — > cr + 7 can be written in the following 
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form: 



M 



M 



e* -U k-p-e* -pU ■ k) , 



(2) 



where e M (/c) and U^p) are the four- vectors describing the vector polarization 
(four- momenta) of the photon and of the V — meson. 

Averaging over the V-meson polarizations and summing over the photon po- 
larizations, one can find for the width of the decay V — > (7 + 7 (still considering 
the a-meson as a stable particle): 

r<V - <n) = .^f* f 1 " ^) 3 - 231 »^ (l - j£) 3 keV, (3) 



where M(m a ) is the mass of the V(cr)-meson. The numerical estimation (3) 
is done for M=0.77 GeV (p- meson). Evidently, at a given value of the gy^- 
coupling constant, the value of T(V — > (77) depends essentially on the value 
taken for the a-meson mass, which can be in a wide range: m a = (0.4 1.2) 
GeV [20]. In Table 1 we report the values of T(p° — > 177) and r(p° — > (77 — > 
7T 7r°7), calculated for the values of the g pcr7 -constant and m c from Refs. [2] 
and [13]. Such estimations for r(p° — > (77 — > 7r°7r°7) are almost two orders of 
magnitude larger than the value given by the experiment [22]: 

T(p° -> (77 -> 7T°A) = (2.85±};f ± 0.6) keV. (4) 



Note that the predictions of QCD sum rules [18] are also in serious disagree- 
ment with the direct experimental estimation (4). 



[GeV] 




iV ^7) 

[keV] 


r(p° 7T 7r°7) 

[keV] 


Ref. 


0.5 


3 


379 


126 


[2] 


0.6 


5-6 


308-444 


103-148 


[13] 



Table 1 



Radiative widths for different masses of a- meson, m a and coupling constant g po - 7 . 

This crude preliminary estimation leads to a conclusion which is inconsistent 
with the suggested interpretation of the HERMES effect [13] and the descrip- 
tion of threshold behavior of the cross section for 7 + p — > p + p° [2]. Therefore 
let us estimate more precisely the g paJ coupling constant on the basis of the 
decay p° — > (77 — > 7r°7r°7, removing the hypothesis taken above, i.e. taking 
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into account the finite value of the a-width, and the possibility that in some 
cases it is possible to have m CT > M. 

Considering the sequence of decays p° — > a + 7 — > 7r° + 7r° + 7, one can find, for 
the effective n° + 7r°-mass distribution, the following formula (see Appendix): 

dr e 2 g 2 pal ( w 2 \ 3 p w m a Ya 

dw 2 288tt 2 ^ M 2 ) (3 a (w 2 -ml) 2 + Tlm 2 a : U 



where T a is the total width for the a-meson, and: 



1-^ 

1 9 1 yv 



1 - 



4m 2 



For the estimation of the coupling constant g 1 on the basis of r(p° 
7r°7r°7) the following formula can be used: 



(77 



2 r(p° - > a 7 -> 7T 7r°7)727r 

9 pay 



Ma 



nr 2 f(r 1 ,r 2 ) 



(6) 



with 

f ^ )= l d \x% 2 lr 2 r 2 

Using the experimental information for the branching ratio for the radiative 
decay p° — > 0-7 — > 7r° + n° + 7 [22], with the help of Eq. (6), one can deduce 
the constant g paJ as a function of the cr-meson parameters, its mass and its 
width (Fig. 1). One can see a strong dependence of the value of g pay on the 
two coordinates r\ and r 2 , which vary in the wide interval suggested by [20]. 
Moreover, in the whole region of parameters 77 and r 2 one can see that g 1 < 
2.5; for example, for m a = 0.5 GeV and T a = 0.6 GeV, we find g 2 = 0.196. 



1 — a/x 



M 



r 2 



M' 



4m 2 
M2" 



3 Estimation of g VaJ couplings from a — > 27 



Another source for an independent estimation of the 5f pcr7 coupling constant, 
which confirms our previous finding, is the radiative decay a — > 27. In principle 
it is possible to estimate the width of this decay analyzing the a-contribution 
in the amplitude of the elastic scattering of photons by nucleons, in the frame- 
work of the t— channel exchange, Fig. 2a, or in the two-photon production of 
a 7r + 7r~ or 7r°7r°-pair in 77-collisions, Fig. 2b, in the corresponding region of 
the 7T7r-effective mass [24]. 
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0.6 0-8 



Fig. 1. Two-dimensional plot of the coupling constant g pai as a function of r\ = -j^- 
and r 2 = ^. 




Fig. 2. Feynman diagrams for er-meson contribution to nucleon Compton scattering 
(a) and to 7r + 7r~-production in 77-collisions (b). 

Following the VDM approach, (Fig. 3), the matrix element for a — > 27 can be 
written in terms of the gy al coupling constants, as follows: 

M(a -> 27) = e 2 (g p g pal + &,0 W(ry )ei ■ e 2 m 2 a /M, (7) 



where e.\ and e 2 are the three- vectors of polarization for the two produced pho- 
tons. The coupling constant gy determines the width of the V^-meson leptonic 
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7 r(k 2 ) r(k 2 ) 



Fig. 3. Feynman diagrams for a — > 27-decay in the VDM. 

4:71 

decay V — > e + e~ (through one-photon exchange): r(V — > e + e~) = a 2 gyM— . 

o 

Note that this formula holds for zero lepton mass. The existing data about 
p{uj) — > e + e~ decays [20] allow to estimate g p and g w with good accuracy: 
^ ~ 0.04, <£~0.(%J. 

As both couplings g pal and g wo - 7 enter in Eq. (7), we assume, for simplicity, 
the validity of the SU(3) relation: g^/'Wy — 3. The width T(a — > 27) can 
be written as: 

r ( . 2o0 = (l + 13.7^ ( r ^ v ) 3 keV. (8) 



Taking g p(T7 ~ 3 [2], one finds 

r(o- -> 2 7 ) = 123 keV, for m a = 1 GeV. (9) 



The same or even larger numbers for T(a — > 27) follow from the QCD-sum rule 
estimation for g pal [18] or from the value deduced from the analysis of p° — > 
7r + 7r - 7-decays [11]. These numbers contradict the corresponding experimental 
estimations: 

r(<7 -> 2 7 ) = 10.6 keV [20], r(cr -> 2 7 ) = (3.8 ± 1.5) keV [24]. (10) 



And, again, the discrepancy is large, up to one order magnitude. The situation 
is even worse, taking the value of g pal from Ref. [13], which leads to 

r(<7 -> 2 7 ) = (340 -T- 490) keV. (11) 
On the other hand, taking the experimental data about the a — > 27-decay, 



7 



(10), one can deduce directly a value 



) 



3 



0-28 < gl 



pal U GeV 



< 0.78, 



(12) 



in evident contradiction with the previous large g pal values [2,11,13,18], but 
in good agreement with the present estimation, based on the new data about 
the decay p -> 7r°7r°7 [22]. 



4 Discussion and conclusions 

Let us discuss the possible consequences of a smaller value of the g pal coupling 
constant; more precisely, let us consider the case g 2 < 1. Here this range is 
derived from two independent sources: the radiative decay of the p°-meson, 
p° — > o"7 — > 7T° + 7T° + 7, on one side, and the radiative decay of the cr-meson: 
a — > 27, from another side. Both these decays can be described in a relative 
simple and transparent theoretical framework: the cr-dominance for the first 
one and the VDM approach for the second one. Note that the VDM model 
gives a good description of the different numerous radiative decays involving 
vector and pseudoscalar mesons. The coupling constants g p and g u (for the 7 — > 
^-transition), which enter in our consideration of the decay o — > 27, are well 
known from the experimental data about the decays V — > £ + £~ , which have 
quite good accuracy. Therefore, the main uncertainty in the determination of 
g pal on the basis of the existing data on a — > 27, derives from the relatively 
large interval for Yio — > 27). The parameters of the cr-meson also do not affect 
very much this estimation: there is no dependence on the a-width, and only 
a cubic dependence on the a— meson mass, so that: 



where (<?L. 7 )o is the value of the considered coupling constant for m a —l GeV. 
From this scaling law it follows the g pal constant gets smaller for higher a- 
masses. 

The estimation of g 2 pai on the basis of the decay p° — > a + 7 — > 7r° + 7r° + 7 
depends essentially on the mass and width of the cr-meson, only, without 
additional unknown coupling constants. We can build a two-dimensional rep- 
resentation on g pal , with non trivial dependence on m a and Y a . It turns out 
that, for the 'standard' value m a =0.5 GeV, in the interval 0.6 < T a < 1 GeV, 
one has the following limits: 




(13) 



0.19 < g* < 0.30, 



(14) 
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definitely smaller than the values used in [2] and [13,18]. On the other hand, 
the interval (14) is in agreement with the interval (12). We must stress that 
the value = 9 in the consideration of the process 7 + p — > p + p° [2], 
has been found for m CT =0.5 GeV, which has been chosen quite arbitrarily. 
The a-exchange amplitude and the corresponding estimation of g pal depend 
essentially on m a . For example, in the near-threshold conditions for 7 + p — > 
p + p°, another scaling law holds for g pal : 

, 2^ x 2 = const, (15) 
(m* + a) 2 



where a = M 2 m N / \m N + M) ~ 0.32 GeV 2 (m N is the nucleon mass). This 
relation shows a large correlation between g pal deduced from the data about 
7+p — > p+p° and the cx-mass, in such a way that larger value of m a correspond 
to larger value of the coupling constant. 

Moreover, the differential cross section for 7 + p — > p + p°, calculated in frame- 
work of a-exchange, contains a strong dependence on the cr-mass, through the 
phenomenological form factor, which has to be introduced here to improve the 
t— behavior of the differential cross section dai^p — > pp°)/dt at large \t\. To 
deduce a definite value of the g pal coupling constant from a fit of the differ- 
ential cross section data, it is also necessary to assume a definite value for the 
go-AW-coupling. Typically this value is determined from the TVTV-potential. 

The estimation of g pal from the data about p° — > n + + 7r _ +7, which suggests 
the large value: g 2 > 10, can not be considered very precise, as the possible 
a-exchange is hidden by the large contribution of bremsstrahlung (photon 
radiation by charged pions). Therefore the coupling g pai derived by this decay, 
results in a very large width for T(p° — > 7r° + 7T° + 7), in strong contradiction 
with the experiment and and with the other theoretical estimations of this 
width. 

The HERMES-effect can not be considered a reliable source of information 
about the g po - 7 -coupling constant. It is more correct to say, that, in order to 
explain this effect in the framework of the model [13], one needs large values, 
which are difficult to justify on the basis of what is known from the radiative 
decays V — > nnj. For example, the large value of g wul was justified in [13] by 
taking the lower limit for uo — > 7r + 7r~7 from [20], instead of a much smaller 
and more precise value from u — > 7r°7r°7. Following selection rules, the relation 
Y(uo — > 7r + 7r~7) = 2Y(uo — > 7r°7r°7) [23] is correct independently on the decay 
mechanism. One deduces, for gw 7 , a value which is two orders of magnitude 
lower then in Ref. [13]. 

If the small values for the g po - 7 -coupling constant, as we suggest in the present 
paper are correct, the interpretation of the following different electromagnetic 
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processes has to be revised: 

• The explanation of the near-threshold cross section for 7 + p — > p + p° 
in frame of the a-model can not work, because the a-contribution has to 
be small. The a-contribution, taking g 2 < 1, turns out to be one order of 
magnitude smaller than necessary for the explanation of the existing data. In 
principle, one can recover by increasing correspondingly the go-ATjv-coupling 
constant, but this would seriously modify the iViV-potential. Moreover, such 
increase would change the scale of the a-contribution to the differential cross 
section ot^ + p^p + u, in the near threshold region, making the a and 
7r-exchanges of the same order. So other contributions, such as Pomeron 
and (or) / 2 -exchanges, which have been extrapolated up to threshold, have 
to be taken into account. This allows to decrease g pal up to unity [5]. 

• The explanation of the HERMES effect, in terms of coherent contribution 
of mesons, does not hold anymore. 

• The predictions of these constants in the framework of the QCD-sum rules 
do not hold. 
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6 Appendix 

In this Appendix we discuss the effect of the a-width on the decay p — > 

7T + 7T + 7. 

The matrix element, corresponding to this diagram (Fig. 4) can be written as: 
M = M(Va*-f)^— — -m a g U7Tn , (16) 

where A4(Vcr*'y) is the matrix element for the decay V — > a* + 7, with pro- 
duction of a virtual a-meson, which mass w (different from the mass of the 
a-meson, rria), coincides here with the effective mass of the produced it + -n~- 
system, m a — > m CT — iT a /2, T a is the total width of the a-meson. 

Taking the expression for M.{Va*^), Eq.(2), (where the gy CT * 7 coupling con- 
stant does not depend, following our assumption, on the virtuality of the 
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r (k) 




V(q) n-(k 2 ) 



n + (k,)\ 



Fig. 4. Feynman diagrams for V — > 7r + +7r +7 decay through cr-meson production. 

cr-meson), one can find, after summing over the photon polarizations and 
summing over the polarization of the V— meson: 

IMV^D^^-Mlgdl, (17) 



We are interested here in the cir/eiw 2 -distribution (integrated over the effective 
7r + 7r~-mass) for the decay V° — > n + + 7r~ + 7 (we are considering here 7r + 7r~- 
production to have different final pions): 

[A4| 2 , 4 oPfc [ d?ki d 3 k 2 

dT = w (27r) (2^; y (2^(2^ ^ 

As |.M| 2 depends only on the variable w 2 , it can be taken outside the integra- 
tion, so that: 

_ \M\ 2 d*h 

dT ~ 2VME; 1 ' (18) 



where X is the result of an invariant integration (which has to be done in the 
CMS of the produced 7r + 7r~-system): 



1= I ——6( q -h-k 2 ) = 2n-, E* = -,f = ^ T -ml (19) 



P / 

So X = 2nf3 w , where /3 W = — = \/l — tf/w 2 is the velocity of the pion 

produced in the decay of the virtual a- meson, a* — > 7r + + 7r~, with mass w. It 
is the energy of the photon, produced in the decay V — > 7r + + 7r~ + 7, which 
determines the effective mass of the produced n + + 7r~-system, through the 
following relation: 

M 2 = w 2 + 2ME 1 (20) 
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which holds in the rest system of the decaying V— meson. Therefore we can 
write: 

d 3 k ( w 2 \ 

— -> 4irE.dE, = -ir\l- —J dw 2 . (21) 



Substituting Eqs. (19) and (21) in (17) one can find: 



d£ \M\ 2 ( w 2 \ ml 
dw 2% 3 M I M 2 / V w 2 



9 / 9 \ 3 / 9 9 

w \ I mt o m 



^(1-^1 \i l -^9i^- —o ZZ2T2 ( 22 ) 



37r 3 2 9 ^ M 2 y V w 2 ^ w ^ (w 2 - ml ) 



The coupling constant 5W+7T- can be related to the total width IV of the a- 
meson, assuming that a — > 27T is its main decay: r CT = T(a — > 7r + 7r~) + r(cr — > 
7T°7r ). In terms of the coupling constant g = g aw + w - = g a ^^ (on the basis of 
isotopic invariance), one can write: 



4- x \M(a — > 7T+7T") 2 fl 2 „ 

r(<7 - TT+vr-) = — ] -q = ^-m af 3 a , (23) 



87im a 16n 



where f3 a = yjl — Am 2 / m 2 is the velocity of the pion produced in the decay 
of the real a- meson, with mass m a . Taking into account the relation: T(a — > 
7T + 7r _ ) = 2T(a — > 7r°7T ), due to the identity of the two neutral pions in the 
decay a — > 7r°7r°, one can write: 

TV = \v{a -> 7T+7T-) = |-^m<A, (24) 

2 o2 7T 



Substituting Eq. (24) in (22) one obtains: 



3 

dr e 2 g Va * y ( w 2 \ p w m 7T Y a 



dw 2 U4ir 2 V M V (w 2 - m 2 ) 2 + T^m 2 



So, for the full width of the decay V — > 7r + + 7r +7, the following formula 
holds: 

M 2 m 

/• dr 

r(v -> 7T + 7r- 7 ) = / -f^^ 2 - (26) 

j aw 2 

4m 2 
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For the decay V — > n° + it + 7, Eq. (26) holds, divided by a coefficient 2 - 
due to the identity of pions in the final state. 

Let us test the validity of these formulas, considering the limit of zero width, 
F a — > 0, i.e. considering the a-meson as a stable particle. We use the following 
symbolic relation: 



r CT ->o 



(w 2 — ml) — iY a m a 
1 

{w 2 — m 2 ) + iT a m a 



itr 



ml) + iT a m a (w 2 - ml) - iY a m a 



w 2 — m 2 



+ i7i5(w 2 — m 2 ) 



7T 



S(w 2 — m 2 a ) 



This 5-function allows one to integrated easily Eq. (26), getting the following 
result: 



nV - oi - vr+vr-7) = -g^M 1 - ^ 



(27) 



To find the radiative width for the radiative decay V — > (77, taking into 
account both possibilities, a — > 7r + + 7r _ and cr — > 7r° + 7r°, we have to introduce 
a factor 3/2, so that: 

rC-<nr) = ^^(i-^)' 

in agreement with the direct calculation, see Eq. (3). 
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